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' Abstract 

Partially quenched Quantum Chromodynamics with Wilson fermions on a lattice is 
considered in the framework of chiral perturbation theory in the mesonic sector. Two 
degenerate quark flavours are associated with a chirally twisted mass term. The pion 
masses and decay constants are calculated in next-to-leading order including terms linear 
in the lattice spacing a. 

PACS numbers: ll.15.Ha, 12.38. Gc, 12.39.Fe 
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A notorious problem for numerical simulations of Quantum Chromodynamics with dynam- 
ical quarks is to reach the region of realistic small quark masses. To deal with this problem 
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different methods are being used. On the theoretical side, chiral perturbation theory sup- 
plies us with formulae which allow to extrapolate pion masses and other observables from 
medium to small quark masses EJ E]. The low energy parameters of chiral perturbation 
theory, the Gasser-Leutwyler coefficients, can in turn be determined by numerical simulations 
of lattice QCD, see |H for a review. The effects of the finite lattice spacing a are taken into 
account in chiral perturbation theory in form of an expansion in powers of a El 13 IHl El ■ In 
[inillllllHIIHllIlllin] numerical results from lattice QCD are compared with chiral perturbation 
theory in next-to-leading order. 

On the algorithmic side partially quenched QCD is an approach to the regime of small quark 
masses. In a numerical simulation of partially quenched QCD the Monte Carlo updates are 
being made with sea-quarks, which have large enough masses ms in order to allow a tolerable 
simulation speed. On the other hand, quark propagators and related observables are evaluated 
with smaller valence-quark masses my. Chiral perturbation theory has been adopted to the 
case of partially quenched QCD in [1?^ IT7| . 

Recently it has been advocated to employ a chirally twisted quark mass matrix for Wilson 
fermions [iHlIini in order to improve the efficiency of QCD simulations [201 El]. At present, this 
approach appears to be very attractive for reducing lattice effects in numerical simulations. First 
numerical results using this approach are very promising [221 • The corresponding modifications 
of chiral perturbation theory for lattice QCD have been worked out to order a in [23] • 

Combining these approaches, which is a logical next step, it appears attractive to simulate 
QCD with chirally twisted quark masses in a partially quenched manner. For the theoretical 
analysis of the data it is desirable to extend the results of [2H] to the partially quenched case. 
This is the object of this letter. 

We consider lattice QCD with Nf = 2 flavours of sea quarks and the same number of 
valence quarks. For simplicity we restrict ourselves to the case of degenerate quark masses 
(■rriu = ma = m). The quark mass matrix contains two chiral twist angles. Chiral perturbation 
theory is applied to the mesonic sector. The pion masses and decay constants are calculated in 
chiral perturbation theory in next-to-leading (one-loop) order, including lattice terms linear in 
the lattice spacing a. 

Recent Monte Carlo simulations of the partially quenched model [HI [11] show that the 
one-loop contributions of chiral perturbation theory (for vanishing twist) can be identified in 
the numerical data and used to estimate Gasser-Leutwyler coefficients. In particular, it turned 
out that partially quenched chiral perturbation theory for lattice QCD is already applicable for 
quark masses below \ms. 

A theoretical description of partially quenched QCD can be obtained through the intro- 
duction of ghost quarks [21]. For each valence quark a corresponding bosonic ghost quark is 
added to the model. The functional integral over the ghost quark fields then cancels the fermion 
determinant of the valence quarks and only the sea quark determinant remains in the measure. 
In our case there are 2 flavours of valence quarks, sea quarks and ghost quarks, each. The 
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quark mass matrix is 

M = diag(mv, my, m^, ms, mv, mv). (1) 
A chirally twisted mass matrix, depending on two twist angles can be introduced in the form 



M{ljv,ujs) 



where 



(2) 



(3) 



In chiral perturbation theory the dynamics of the pseudo-Goldstone fields is described by 
a low energy effective Lagrangian. For partially quenched QCD the pseudo-Goldstone fields, 
which we shall generally call pions, are parameterized by a graded matrix of the form 
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The 4x4 matrix A and the 2x2 matrix D contain commuting numbers, whereas the 4x2 
matrix B and the 2x4 matrix C have anticommuting entries. U is an element of the supergroup 
SU(4|2) and can be represented as 



U{x) 



exp 



(5) 



The commuting elements of the graded matrix $ represent the pseudo-Goldstone bosons made 
from a quark and an anti-quark with equal statistics, and the anticommuting elements of $ 
represent pseudo-Goldstone fermions which are built from one fermionic quark and one bosonic 
quark. The supertrace of $ has to vanish, 



^r $ = 



(6) 



where one defines 



sTr 



TvA - TvD . 



(7) 



A B 
C D 

This condition introduces technical complications in loop calculations. One way to deal with 
them [in] is to drop the condition sTr$ = 0, adding a mass term proportional to (^r$)^ and 
sending its mass to infinity in the end. Another way, which we followed in our calculations, is 
to choose a basis of 35 super-traceless generators Tj of SU(4|2) and to expand 



$ = ^7r,T,. 



(8) 



The chiral effective Lagrangian to leading order, including lattice artifacts up to first order 
in the lattice spacing a is given by [7] 



rp2 



(9) 



The mass term contains the matrix 



and the lattice term is parameterized by 



X 



2BqM 



p = 2Woa 1 



(10) 



(11) 



with two constants Bq and Wq. 

Similar to the unquenched case, the chiral twist can be incorporated in chiral perturbation 
theory by transforming the lattice term into 
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The effective Lagrangian in next-to-leading order (with respect to masses and momenta in the 
usual Weinberg power counting scheme) is then given by [Z| 



sTr [d^Ud^U^) - ^ sTr (xU^ + Ux^) - ^ sTr {p{u)U^ + U p{ 



(13) 



sTr (d^Ud^U^ 



d„Ud„U^ 



L2 sTr (d^Ud.U^) sTr [d^Ud.U^ 
- U sTr (dMdJJ^) sTr (xU^ + Ux 



+Wi ^r (d^Ud^U^) ^r [p{uj)U^ + U p{uj)^) + L5 sTr (d^Ud^U^ [xU^ + Ux^ 

sTr (xU^ + Ux' 



sTr (d^Ud^U^ [p{^)U^ + Up{uj\ 
-Wq sTr (xU^ + Ux^) sTr [p{uj)U^ + Up^uj)^ 
-Wj sTr (x?/"^ - Ux^) sTr (p(cj)f/"^ - Up{ 
-Ws sTr (xt/V(c^)f/^ + Up{uj)^Ux^) + ^(a^ 



sTr (xU^ - Ux^ 
LssTv(xUhU^ + Ux^Ux^ 



In order to calculate the pion masses and other physical quantities, the Lagrangian has to be 
expanded in terms of the fields contained in In the untwisted case the expansion is around 
the minimum (more precisely, the saddle point) at vanishing $. The chiral twist, however, 
introduces a shift of the expansion point, as was already observed in the unquenched case 
[28] . We have calculated this shift in next-to-leading order including terms linear in the lattice 
spacing a, and expanded the effective action around the shifted minimum. 

From the resulting expression the propagators and pion decay constants can be calculated 
in a one-loop calculation. The results are known for the case of vanishing twist f^. Therefore 
one has to account for the effects of the twist angles now. The lattice term ^i{p{u)W + U pi^uY) 
yields the quadratic contribution sTr((p(ci;) + p(ci;)^)$^), which contains the combination 



p{uj) + p{iuy = AWoa{cosuJv'l +003^0151 +costi;yl 



G\ 
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We see that in this term the twist amounts to a multiplication of the untwisted lattice terms 
by factors cosuv and cosus, respectively. We have checked all terms proportional to the 
coefficients Wi and found that for them the same rule holds. 



In the one-loop calculation the propagators and vertices from the leading order effective 
Lagrangian enter. We considered these terms and found that for the propagator and 4-point 
vertices the twist produces cos c<j-factors as above. The 3-point vertices are modified by sinu- 
factors. The relevant Feynman diagrams, however, contain two of these vertices and thus 
contribute to order only. 

To summarize, at order a the chiral twist amounts to adding factors cosuy or COSC1J5 to the 
untwisted quantities. This rule will not apply to order a^, where new vertices appear. Results 
for masses and decay constants for partially quenched QCD with Nf = 3 flavours have been 
given by Rupak and Shoresh [Zj. The formulae for a general number Ns of sea quarks and 
Ny = 2 valence quarks are presented in [13]. They give masses and decay constants for flavour- 
charged pious. From these we obtain the results for the twisted case under consideration using 
the procedure explained above. For the the pion masses we find 



Xs + pcosus 



^ss = Xs + P cos Us 

^0 ^0 ^0 
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+^^2(^6 - W4)xspcosuJs + t;2(2^8 - - L5)xspcosiJs 

^0 -f^n 



L4)xsP COS Us 



(15) 



IT^VV = XV + PCOSUJv 

1 



327r2K2 



iXV + pCOSUJv)<XV - Xs + P cos Uy - p cos LUs 



Xv + pcosuv 

A2 



+ 2xv -Xs + ^pcosujv - pcosus) log 



16 8 16 

+ ^^2(2^6 - Li)xsXv + -^C^Ls - L5)xl + ^^2(^6 - Li)xspcosu;v 
^0 ^0 ^0 



+ -^{W(i - Wi)xvpcosuJs + ^^2(2^^8 -W5- L5)xvpcosuv 



(16) 



^vs = -{xv + Xs + P cos ujy + p cos Us) 



iXv + Xs + pcosuv + p COS Us) ixv + P cos uv) log 



Xv + pcosa;y 

A2 



+-^{2Lq - U)xs{xv + Xs) + -E^{2Ls - L^){xv + XsY 

8 8 
+^(W^6 - L4)xsipcosuv + pcosus) + -^{We - W4){xv 

2 

+^{2W8 -W5- L5)(xv + Xs)(pcosa;y + pcosa;5) 



Xs)p cos Us 



(17) 



and for the decay constants 
Fss _ ^ _ 



{Xs + P cos Us) log 



Xs + P cos Us 
A2 



-Tn -Tn -Tn 



(18) 



vv 



XV + XS + P cos ujv + p cos u;^ log 



+ Xs + pcosujv + pcosujs 



2A2 



+ -=;^^5Xv + t;2^4X5 + T;2^5Pcos^^ + T;2^4pcosu;s 

-Tn -Tn -Tn -Tfl 



(19) 



where 



"^{xs + pcosujs) log 



X5 + pcosa;5 
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+ (xy + X5 + P cos ujv + p cos cus) log 
1 



A2 

XV + Xs + P cos CJy + p COS LUs 



2A2 



"^12871^^2 - X5 + pcosu;y - pcosus - 1^X5 + pcosu;5j log 

2 8 2 ^ 

+T;2^5(Xy + Xs) + -B2L4XS + -^W^ip cos ujv + P cos ujs) + 
^0 ^0 ^0 



Xv + pCOSUy 
Xs + pCOSiUs 



^Wipcosus, (20) 



Xv = 2Bomv , Xs = ^Boms , p = 2Woa , (21) 

and A is the renormalization scale. It should be noted that the axial vector current used for the 
definition of the pion decay constant is the Noether current corresponding to an axial rotation, 
which differs from other commonly used axial vector currents by 0{a) terms. 

For maximal twist, uy = ujs = 7r/2, the lattice artifacts vanish to order a, as has been shown 
for lattice QCD in general by Frezzotti and Rossi [213 EI] and checked in chiral perturbation in 
[23] . Terms of order could be incorporated along the lines of [H]. Such a calculation, which 
is technically rather involved, is in progress. 

The expressions for the pion masses and decay constants in partially quenched lattice QCD 
with a chiral twist can be used in the analysis of numerical results from Monte Carlo calculations 
and will aid the extrapolation to small quark masses. 

We thank I. Montvay, F. Farchioni and R. Frezzotti for discussions. 
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